The X— ray Source at the Center of 
the Cassiopeia A Supernova Remnant 



S. Mereghetti 

Istituto di Fisica Cosmica G.Occhialini - CNR 
v.Bassini 15, 1-20133 Milano, Italy 
sandro @ifctr. mi. cnr. it 

A. Tiengo 

XMM-Newton Science Operation Center 
VILSPA ESA, Apartado 50727, 28080 Madrid, Spain 
tiengo @xmm. vilspa. esa.es 

G.L. Israel 

Osservatorio Astronomico di Roma 
v.Frascati 33, I-OOO4O Monteporzio C atone, Italy 
gianluca@mporzio. astro.it 

ABSTRACT 

We present the first results of an XMM-Newton observation of the central 
X-ray source in the Gas A supernova remnant. The spectrum can be fit equally 
well with an absorbed steep power law {uph ~3, N^^ ~1.5xl0^^ cm~^) or with a 
bremsstrahlung with temperature kT~2.4 keV and N// ~10^^ cm^^. A blackbody 
model (kT^B ~0.7 keV) gives a slightly worse fit and requires a column density 
Njif ~5xl0^^ cm~^ and an emitting area with radius ~0.3 km (for d=3.4 kpc). 

A search for pulsations for periods longer than 150 ms gave negative results. 
The 3(7 upper limits on the pulsed fraction are ~13% for P>0.3 s and ~7% for 
P>3 s. 

The overall properties of the central X-ray source in Gas A are difficult to 
explain in terms of a rapidly spinning neutron star with a canonical magnetic 
field of ~10^^ G, and are more similar to those of slowly rotating neutron stars 
such as the Anomalous X-ray Pulsars. 



Subject headings: Stars: neutron; X-ray: stars; Supernovae: individual (Gas A) 
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1. Introduction 

The long-sought compact remnant born in the Cas A supernova explosion was discovered 
in the first light image obtained with the Chandra X-ray satellite (Tananbaum 1999). For 
the following reasons there is little doubt that CXO J232327.8+584842 (CXO for brevity) 
is indeed associated to Cas A , the youngest (~320 yrs) supernova remnant observed in the 
Galaxy (Ashworth 1980, Fesen et al. 1987). 

First, the X-ray fiux of ~10~^^ erg cm~^ s~^ (Pavlov et al. 2000, Chakrabarty et al. 
2001) and the strong upper limits on any optical/IR counterpart in the small (1" radius) 
Chandra error circle (R>26.3, K>21.2, Ryan et al. (2001), Kaplan et al. (2001)), imply for 
CXO an X-ray to optical flux ratio greater than ~800. This immediately excludes any kind 
of non-compact galactic source as well as accreting binary systems, unless the companion 
star is a very low mass object. On the other hand the X-ray to optical flux ratio of CXO 
is similar to that of radio pulsars and other neutron star candidates (see, e.g.. Shearer & 
Golden 2001, Treves et al. 2000). 

Second, CXO is located at an angular distance <7" from the site of the Cas A supernova 
explosion, that has been recently re-determined with good accuracy by a backward projection 
of the expanding optical filaments (Thorstensen et al. 2001). For an age of 320 yrs and 
a distance of 3.4 kpc (Reed et al. 1995), the implied transverse speed of the neutron star 
would be ~340 km s~^, well within the range of the observed space velocities of radio pulsars 
(Cordes & Chernoff 1998). 

Finally, the possibility of a background extra-galactic source has been dismissed (Pavlov 
et al. 2000) on the basis of (a) the very small probability of finding an AGN as bright as 
CXO within an area of only ~0.05 arcmin^, and (b) the steep power law X-ray spectrum of 
CXO. 

Although the association with Cas A is well established, the exact nature of CXO 
remains puzzling. Deep radio observations (McLaughlin et al. 2001) were unable to detect 
any source, either pulsed (upper limit of 1.3 mJy at 1.4 GHz) or unpulsed (u.l. of 40 and 
6 mJy at 1.3 and 4.4 GHz, respectively). A radio/X-ray synchrotron nebula could indicate 
the presence of a rapidly rotating, magnetized neutron star, even when beam orientation 
effects prevent its direct observation as a radio pulsar. However, no radio or X-ray (Murray 
et al. 2001) nebula has been found around CXO. 

The X-ray spectrum of CXO is not well constrained. Equally good fits to the Chan- 
dra data were obtained with simple power law, blackbody and bremsstrahlung models 
(Pavlov et al. 2000, Chakrabarty et al. 2001, Murray et al. 2001). Power law fits give 
a steep photon index Q;p/j~3-4 and a column density N/^~(1.5-2.3) xlO^^ cm~^, larger than 
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the estimated value of the interstellar absorption in the direction of Cas A (Keohane et al. 
1996). Young neutron stars emitting in the X-ray band by non-thermal magnetospheric 
processes have harder spectra («p^~1.5) and higher luminosity (Lx~10^^-10^^ erg s~^) than 
CXO. The fits with a blackbody model give kT^B^O.S keV and N//~10^^ cm"^, but they 
imply an emitting area that, for any reasonable assumed distance, is far too small to be con- 
sistent with the surface of a neutron star. This discrepancy is not uncommon when fitting 
cooling neutron stars with simple blackbody spectra. Larger emitting areas, more consistent 
with the neutron star dimensions, are generally obtained when more appropriate models for 
the thermal emission from neutron star atmospheres are used (see, e.g., Zavlin et al. 1998, 
1999). However, Pavlov et al. (2000) showed that in the case of Cas A also fits with neutron 
star atmosphere models imply emission from only a small fraction of the star surface. 

A re-analysis of archival Einstein and ROSAT observations showed no evidence for long 
term variabihty (Pavlov et al. 2000). Search for X-ray pulsations have been carried out 
with Chandra by Chakrabarty et al. (2001) with negative results, and later by Murray et 
al. (2001) who reported a possible periodicity at 12 ms, but with a very low statistical 
significance. If this period were confirmed, it would be difficult to explain the lack of a 
strong pulsar wind nebula around CXO, unless its magnetic field is rather low (<5xl0^° G). 

Here we report the results of an observation of CXO performed with the XMM-Newton satel- 
lite on July 27*'*, 2000 as part of the Performance Verification program. 

2. Observations 

The data reported here were obtained with the European Photon Imaging Camera 
(EPIC) instrument on the XMM-Newton satellite. EPIC consists of two MOS CCD detectors 
(Turner et al. 2001) and a PN CCD instrument (Striider et al. 2001), giving a total collecting 
area > 2500 cm^ at 1.5 keV. The mirror system offers an on-axis FWHM angular resolution 
of 4-5" and a field of view of 30' diameter. 

The XMM-Newton observation of Cas A started at 6:45 UT of July 27*'*, 2000 and 
lasted about 11 hours. The telemetry saturation due to the high overall counting rate from 
the supernova remnant caused a non-continuous operation of the CCDs, yielding effective 
exposures of 9219 s in the PN, and 23573 s and 23594 s in the MOSl and M0S2 detectors, 
respectively. The satellite was pointed at the center of the remnant (RA: 23'* 23"* 25*, Dec: 
58° 48' 20" (J2000)), resulting in the detection of CXO very close to the center of the EPIC 
field of view (off-axis angle of 0.52'). During this observation the thin filter was used. The 
MOS cameras were operated in Large Window mode to limit the source photon pile-up. 
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The background corrected count rates of CXO were 0.25±0.02 and 0.13±0.01 counts s ^, 
respectively in the PN and in each MOS. 

The data were analyzed with version 5.2 of the XMM-Newton Science Analysis System. 
The counts for the spectral analysis were extracted in the three EPIC cameras from circles 
with 10" radii centered at the Chandra position of CXO. The corresponding background 
corrected count rates were 0.25±0.02 and 0.13±0.01 counts s~^, respectively in the PN and 
in each MOS. Such a small radius, enclosing only about 50% of the counts from a point source, 
was chosen to avoid contamination from bright nearby filaments due to the emission from 
the supernova remnant (see the Chandra images of Murray et al. (2001) and Chakrabarty et 
al. (2001)). All the spectral fits reported below were obtained with the appropriate response 
matrices that take into account the dimensions of the extraction region. The derived fluxes 
are corrected for the fraction of counts outside the 10" circle. 

CXO is relatively faint compared to the non-uniform, diffuse emission due to Cas A . 
We estimate that the background accounts for ~75-80% of the counts within the adopted 
extraction radius (14050 for the PN and 13320 for the MOSl). Therefore, we devoted 
particular attention to the choice of the region for the extraction of the background spectra. 
Based on the high resolution Chandra images, we tried different regions with a surface 
brightness similar to that at the position of CXO and avoiding bright filaments. We found 
that the resulting CXO spectra are somewhat dependent on the chosen background that 
can lead to the presence of residuals near the energies of the strongest Cas A emission lines. 
Such residuals in the background subtracted spectra were minimized using the background 
from a rectangular (~20xl5 arcsec^) region to the East of CXO, roughly corresponding to 
the eastern half of that used in the Chandra analysis by Chakrabarty et al. (2001). The 
spectra discussed below were obtained with this choice of the background region. 

For the PN spectra we used photons corresponding to single and double events, after 
checking that the individual single and double events spectra gave the same results within 
the uncertainties. The energy channels were grouped to have at least 30 net counts per 
bin, and the fits were performed in the 0.6-10 keV energy range. The results of power-law, 
blackbody and bremsstrahlung fits are summarized in Table 1. 

The same extraction and binning criteria were applied to the two MOS cameras. While 
the results obtained with MOSl agree with the PN ones (see Table 1), the M0S2 fits gave 
different parameters and worse values. We found that this discrepancy is due to the fact 
that at the time of the XMM-Newton observation of Cas A the offset calculation for the 
MOS cameras was still based on an imperfect on-board algorithm (only later in the mission 
better fixed tables were implemented). This problem can cause spurious event detections in 
channels close to the low energy threshold and also produces an offset of the energy scale 
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in particular regions of the CCD. By carefully examining the M0S2 image and comparing 
spectra of different regions, we verified that in this observation such an effect occurred just 
at the position of CXO in the M0S2. Since these effects cannot be easily corrected, we did 
not consider the M0S2 spectra in our further analysis. 

To better constrain the spectral parameters, we fitted jointly the PN and MOSl data 
(Table 1). The power law and bremsstrahlung models give better fits to the data, while 
the fit with a blackbody (Fig.l) yields higher values. For comparison with the results of 
Murray et al. (2001), we present in Table 1 also the values obtained by fixing Ni:/=l.lxlO^^ 
cm~^, and by adding a power law to the blackbody model. For all the models, the observed 
0.1-10 keV flux from the joint PN-MOSl flts is of the order of ~2xl0~^^ erg cm"^ s"^. The 
systematic difference between the PN and MOSl fluxes must be ascribed to the uncertainties 
in the relative calibration of the two detectors. The errors on the luminosities reported in 
Table 1 have been computed taking into account the uncertainties of the spectral parameters 
and absorption (90% c.l. for two interesting parameters). 

For the timing analysis we extracted all the counts with energy in the 1 to 6 keV 
energy range from circles with radius 10" centered on the CCO position. This yielded 
respectively 12276, 12652 and 12511 counts for the PN, MOSl and M0S2 detectors. The 
upper limits reported below take into account that the background, estimated from nearby 
regions, accounts for ~75% of these counts. After correction of the arrival times to the 
Solar System baryccntcr, a Fourier analysis of the background subtracted light curves was 
performed using the method described in Israel & Stella (1996). The frequency resolution 
was 2.46x10"^ Hz. The frame integration times of 0.073 and 0.9 s, respectively in the PN 
and MOS cameras, set the upper limits to the explored frequency ranges. No signiflcant 
periodicity was found. The upper limits on the pulsed fraction, computed according to 
Vaughan et al. (1994), are shown in Fig. 2 as a function of the frequency. The curves refer 
to the 3(7 upper limit on the source pulsed fraction in the 1-6 keV energy range assuming 
a sinusoidal modulation. Above 0.4 Hz only the PN data could be used, while at lower 
frequencies we could sum the power spectra from the three instruments, thus improving the 
sensitivity as shown by the lower curve of Fig. 2. Also separate searches in the soft (1-2 keV) 
and hard (2-6 keV) energy ranges gave negative results. The corresponding upper limits are 
about a factor two higher than those shown in Fig. 2. 

3. Discussion 

Our spectral results can be compared with those derived from the Chandra observations. 
In the power law case, our best flt values for a^h and N// are consistent with the values found 
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by Chakrabarty et al. (2001), but not with the more recent findings of Murray et al. (2001) 
who obtain a softer and more absorbed spectrum. The temperature values that we find 
for the bremsstrahlung and blackbody models are higher than the Chandra values, and 
the corresponding column densities slightly smaller. The blackbody normalization of our 
fits implies an emitting area with radius of 0.32d3.4fcpc km, smaller than that inferred with 
Chandra . To check whether the discrepancy with respect to the Chandra results might be 
explained by the higher background contribution in the EPIC spectral extraction region, we 
also performed an analysis with smaller source extraction regions, down to a radius of 5" 
(always using the appropriate response matrices). The derived best fit parameters did not 
vary significantly. 

The above mentioned difficulties in interpreting the X-ray emission from CXO in terms 
of a young active pulsar are therefore confirmed by our XMM-Newton results. An indepen- 
dent knowledge of the expected column density could help to discriminate among the various 
spectral models. Values in the range (0.9-1.3) x 10^^ cm~^ have been derived by Keohane et 
al. (1996), with significant variations across the remnant and with relatively smaller absorp- 
tion at the position of CXO. On this basis, the thermal bremsstrahlung fit is the one giving 
the most consistent Nh value. The low absorption, as well as the higher value, suggests 
that the simple blackbody fit is not adequate. The trend of the residuals shown in the lower 
panel of Fig. 1 (see also Fig. 3 of Murray et al. (2001)) probably indicates that atmosphere 
models more complex than simple blackbody spectra are required. 

The soft spectrum and high X-ray to optical flux ratio make CXO similar to the Anoma- 
lous X-ray Pulsars (AXPs, Mereghetti 2001a). These objects have periods in the 6-12 s 
range, steadily increasing on timescales of ~10^-10^ years, and two (possibly three) of them 
are located at the center of shell-like supernova remnants. Similar properties (except for the 
periodicity) are shared by the X-ray sources at the center of the supernova remnants Puppis 
A (Petre et al. 1996, Zavhn et al. 1999), G290.5+10.0 (Mereghetti et al. 1996, Zavhn et al. 
1998), RCW 103 (Gotthelf et al. 1999), and G266.1-1.2 (Slane et al. 2001, Mereghetti 2001b, 
Pavlov et al. 2001a). The bolometric blackbody luminosity and temperature of all these 
sources arc plotted in Fig. 3, where for comparison we have also included two classes of older 
neutron stars: three middle aged radio pulsars from which thermal emission from internal 
coohng has been observed (PSR B1055-52, PSR B0656-I-14, and Geminga; see, e.g., Becker & 
Triimper 1997) and several nearby isolated neutron stars (Treves et al. 2000, and references 
therein). For all the sources in which also a non-thermal spectral component is present, we 
have plotted in Fig. 3 only the bolometric luminosity of the blackbody component. 

The AXPs have on average a higher luminosity than the central compact objects (CCOs) 
in supernova remnants, but their blackbody temperatures are very similar. Although the 
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majority of AXPs and CCOs have quite stable luminosity, there are two interesting variable 
sources that might be intermediate objects linking these two classes. The candidate AXP 
AX J1845-03 has been observed since 1997 at a flux level corresponding to Lx^lO^^ erg 
s~^ (Torii et al. 1998, Vasisht et al. 2000), at least ten times smaller than at the time of 
its discovery. IE 1614-5055, the CCO in RCW 103, aready known to be variable between 
Fo.5-2fceV~3 and ~50xl0~^^ erg cm~^ s~^ (Gotthelf et al. 1999) and with a possible ~6 
hr periodicity (Garmire et al. 2000a), has recently been detected with Chandra in an even 
higher state (Garmire et al. 2000b). 

Alpar (2001) proposed a unifled interpretation of the classes of objects shown in Fig. 3, 
in which they are all assumed to be neutron stars with ordinary magnetic fields (10^^-10^^ G), 
thus avoiding the need for super-critical magnetic fields required by the " Magnetar" models 
for AXPs (e.g., Thompson & Duncan 1996). According to Alpar (2001), isolated neutron 
stars experience a mass infiow due to the formation of a residual accretion disk from fallback 
material after the supernova explosion. This causes an evolutionary phase in the propeller 
regime, before the neutron star turns on as a radio pulsar or as an AXP, depending on the 
initial conditions in terms of rotational period, magnetic field and mass in the residual disk 
(see also Chatterjee et al. 2000 for a similar model). The dim, nearby neutron stars in the 
lower left corner of Fig. 3 are interpreted as examples of sources in the propeller regime, in 
which the X-ray emission is from cooling powered by internal friction. The CCOs would be 
the propeller objects with the highest mass inflow, surrounded by an optically thick corona 
and might evolve into AXP (Alpar 2001). 

Searches for periodicities in the CCOs extending to low frequencies are therefore rele- 
vant, since some of these objects might have periods similar to those of the AXP. In this 
respect, we note that, while the CCO in G290. 5-1-10.0 has a fast period of 424 ms (Zavhn 
et al. 2000) and might be a "radio-quiet" pulsar similar to Geminga, the fast periodicity 
at 75 ms reported in the Puppis A CCO (Pavlov et al. 1999) has not been conflrmed by 
more accurate Chandra observations (Pavlov et al. 2001b). For frequencies v <6 Hz, our 
limits on the pulsed fraction of CXO are smaller than those previously obtained with Chan- 
dra (Chakrabarty et al. 2001, Murray et al. 2001). Only one of the AXPs has a pulsed 
fraction as small as 7% (see, e.g., Israel et al. 2001), our upper limit for long periods. 



4. Conclusions 

We have reported the flrst results of an XMM- Newton observation of CXO, the interest- 
ing X-ray point source at the center of the Cas A supernova remnant. The large collecting 
area of the XMM-Newton instruments has allowed us to derive strong upper limits of ~13% 
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for P>0.3 s and ~7% for P>3 s on the pulsed fraction (Scr). 

The spectral analysis confirmed the very soft spectrum of CXO , but yields, in the 
case of thermal models, higher temperatures than those measured earlier with Chandra . 
The source luminosity does not show evidence of long term variations when compared with 
previous measurements. 

As for other candidate neutron stars in supernova remnants older than Cas A , the 
overall properties of CXO are difficult to explain in terms of a rapidly spinning neutron star 
with a canonical magnetic field of ~10^^ G, and are more similar to those of slowly rotating 
neutron stars such as the AXPs. 

We acknowledge the finantial support of the Italian Space Agency. This work is based 
on observations obtained with XMM-Newton , an ESA science mission with instruments and 
contributions directly funded by ESA Member States and NASA. 
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Table 1. Results of spectral fits 



Model 






kT« 


Flux* 




Xled/dof 




(10^2 cm-2) 




(keV) 


(erg s-icm~2) 


(lO^^erg s-i) 




PN 

Power law 
jjidCKDoay 
Bremsstr. 


1 47+0.34 
-0.27 

n 4fi+0-19 
U-4D_o.i5 

-'-■'-'^-0.18 


9 or+0.31 


U. ( OlCU.UO 
^■"^-0.50 


2.3 10-12 

Z.U iU 

2.2 10-12 


5111°! 

O.OZCU.O 
8-0-1.3 


1.09/30 

i.OO / ou 
1 1 o /on 

l.ld/oO 


MOSl 
Power law 
Blackbody 
Bremsstr. 


1 p-q+0.42 
-'-•^'-'-0.33 


'-'•-'-^-0.39 


0.65±0.07 


1.9 10-12 
1.67 10-12 
1.8 10-12 


00+364 
q 1 +0.5 
Q 9+3.5 


1.04/29 
1.27/29 
1.03/29 


PN+MOSl 

Power law 

Blackbody 

Bremsstr. 

Power law 

Blackbody 

PL+bbody 


J^-'-'<J-0.14 

1.1 (fixed) 
1.1 (fixed) 
1.1 (fixed) 


9 nfj+0.25 
-i.yu_o.23 

2.60±0.11 

9 an +0.33 
■^■"-^ -0.32 


0.70±0.05 

9 qc+0.44 
.i.ou_o.36 

0.60±0.03 

0.62 tr^i 


2.1 10-12 
1.8 10-12 

2.0 10-12 

2.1 10-12 
1.7 10-12 
2.1 10-12 


q 9+0.3 
'^■^-0.2 

24.3+tj 
4.2l°:l 
19.5l22g« 


1.07/61 
1.49/61 
1.11/61 
1.20/62 
2.00/62 
1.09/60 


a 


Errors at 90% confidence level for a single interesting 


parameter 





* Absorbed flux 0.1-10 keV 
0.1-10 keV luminosity (d=3.4 kpc) corrected for absorption. 



- 12 - 



Fig. 1.— EPIC PN and MOSl spectra of CXO fitted with a blackbody model. 

Fig. 2. — Upper limits (Sex) on the CXO pulsed fraction (semiamplitude of a sinusoidal 
modulation after background subtraction). The curve for v <0.4 Hz was obtained from the 
sum of MOSl, M0S2 and PN data, while at higher frequencies only the PN data could be 
used. 

Fig. 3. — "Color-Magnitude" diagram for different classes of neutron stars: Anomalous X- 
ray Pulsars (stars), supernova remnant Central Compact Objects (squares), nearby Isolated 
Neutron Stars (circles) and radio pulsars (diamonds) . For the luminosity of AXP and radio 
pulsars only the thermal spectral component is considered. To allow the comparison between 
the different objects we plotted the temperatures inferred from pure blackbody fits; the use 
of atmosphere models would move all the points to lower temperatures by a factor ~2-3. 
The luminosity error bars reflect the distance uncertainties. The two positions for RCW 103 
and AX J1845-03 indicate the observed variability range. 
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